Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,736,914 T→A 5.7% intergenic (+31/‑240) bamD → / → raiA BamABCDE complex OM biogenesis lipoprotein/cold shock protein associated with 30S ribosomal subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,736,9140TA 5.7% 488.4 / 15.0 140intergenic (+31/‑240)bamD/raiABamABCDE complex OM biogenesis lipoprotein/cold shock protein associated with 30S ribosomal subunit
Reads supporting (aligned to +/- strand):  ref base T (75/57);  new base A (2/6);  total (77/63)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 141 total aligned reads displayed.

ACCCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACAAA  >  NC_000913/2736767‑2737061
                                                                                                                                                    |                                                                                                                                                   
aCCCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTg                                                                                                                                                    >  3:397945/1‑149 (MQ=255)
  ccAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAg                                                                                                                                                  >  1:1074592/1‑149 (MQ=255)
   cAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGc                                                                                                                                                 >  1:1170430/1‑149 (MQ=255)
    aGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCt                                                                                                                                                >  3:1065320/1‑149 (MQ=255)
    aGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCt                                                                                                                                                <  2:1470834/149‑1 (MQ=255)
    aGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCt                                                                                                                                                >  2:1414294/1‑149 (MQ=255)
      gCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGc                                                                                                                                              >  1:1290569/1‑149 (MQ=255)
          cgcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCAAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCTCAAGGGCTGCCGtt                                                                                                                                          <  1:76829/149‑1 (MQ=255)
           gcgTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGttt                                                                                                                                         >  3:625732/1‑149 (MQ=255)
           gcgTGATGAGCTGCGGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACGTAACAGAAACCTGAAACAC‑AAAACGGCAGCTCAAGGGCTGCCGttt                                                                                                                                         <  4:966302/149‑1 (MQ=255)
             gTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGttttt                                                                                                                                       >  4:535136/1‑149 (MQ=255)
              tGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGtttttt                                                                                                                                      <  3:1107205/149‑1 (MQ=255)
              tGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGtttttt                                                                                                                                      <  2:820445/149‑1 (MQ=255)
                  gcgcTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTAtt                                                                                                                                  <  1:1053038/149‑1 (MQ=255)
                    gcTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACACAAAAACGGCAGCTCAAGGGCTGCCGttttgtgtttc                                                                                                                                 <  2:1487947/149‑8 (MQ=38)
                      tGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                              <  4:1008357/149‑1 (MQ=255)
                       gCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTc                                                                                                                             <  1:815405/149‑1 (MQ=255)
                       gCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTc                                                                                                                             <  1:1414294/149‑1 (MQ=255)
                          gCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGt                                                                                                                          >  2:626524/1‑149 (MQ=255)
                          gCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGt                                                                                                                          >  3:919980/1‑149 (MQ=255)
                            tGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTtg                                                                                                                        >  2:1519929/1‑149 (MQ=255)
                             taTGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTtgt                                                                                                                       <  3:136952/148‑1 (MQ=255)
                                ggAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGaa                                                                                                                    >  1:425920/1‑149 (MQ=255)
                                 gAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                              >  4:250430/1‑138 (MQ=255)
                                 gAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                              <  3:250430/138‑1 (MQ=255)
                                      tGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                              <  4:1422138/149‑1 (MQ=255)
                                      tGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                              <  3:1201031/149‑1 (MQ=255)
                                         aTACCGTCAGATGCAGATTAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATCCATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                              >  2:1176411/1‑146 (MQ=255)
                                         aTACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGaa                                                                                                              <  1:1176411/146‑1 (MQ=255)
                                         aTACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCg                                                                                                           >  2:501873/1‑149 (MQ=255)
                                         aTACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCg                                                                                                           >  2:303305/1‑149 (MQ=255)
                                           aCCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAt                                                                                                         <  2:750122/149‑1 (MQ=255)
                                            ccGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAtt                                                                                                        <  4:247323/149‑1 (MQ=255)
                                               tCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGaaa                                                                                                                   <  3:542280/135‑1 (MQ=255)
                                                 agatgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTc                                                                                                   <  4:1291049/149‑1 (MQ=255)
                                                 agatgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTc                                                                                                   <  3:265857/149‑1 (MQ=255)
                                                    tgcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGt                                                                                                                              <  1:736445/119‑1 (MQ=255)
                                                     gcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTa                                                                                               >  1:940888/1‑149 (MQ=255)
                                                     gcagatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTa                                                                                               >  1:64026/1‑149 (MQ=255)
                                                       agatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGCCATCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGCTGTGATACTGAAGCGATTTAGTCGCTATc                                                                                             >  4:548513/1‑149 (MQ=255)
                                                        gatgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCg                                                                                            >  2:743651/1‑149 (MQ=255)
                                                          tgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGAt                                                                                          <  2:940888/149‑1 (MQ=255)
                                                          tgAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGAt                                                                                          <  3:874338/149‑1 (MQ=255)
                                                           gAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGAtc                                                                                         >  4:899995/1‑149 (MQ=255)
                                                                 gcAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCATTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCtcatca                                                                                   >  4:1247297/1‑149 (MQ=255)
                                                                     gCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAata                                                                               <  2:64026/149‑1 (MQ=255)
                                                                     gCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAata                                                                               >  2:940373/1‑149 (MQ=255)
                                                                      cTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAatat                                                                              <  2:425920/149‑1 (MQ=255)
                                                                             aGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTc                                                                                                                                 <  3:309821/91‑1 (MQ=255)
                                                                             aGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTc                                                                                                                                 >  4:309821/1‑91 (MQ=255)
                                                                                aGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCttt                                                                    >  1:137356/1‑149 (MQ=255)
                                                                                 gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑GAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTg                                                                                                                >  2:841785/1‑104 (MQ=255)
                                                                                 gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑GAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTg                                                                                                                <  1:841785/104‑1 (MQ=255)
                                                                                 gCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTg                                                                   <  4:1021147/149‑1 (MQ=255)
                                                                                   gAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTgag                                                                 >  1:506361/1‑149 (MQ=255)
                                                                                    aaaaTCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTgaga                                                                >  2:511175/1‑149 (MQ=255)
                                                                                     aaaTCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATgg                                                                            >  4:554353/1‑136 (MQ=255)
                                                                                     aaaTCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATgg                                                                            <  3:554353/136‑1 (MQ=255)
                                                                                           tcGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCt                                                         >  1:54581/1‑149 (MQ=255)
                                                                                           tcGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCt                                                         >  2:1268932/1‑149 (MQ=255)
                                                                                              ccgcAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                      <  1:911444/149‑1 (MQ=255)
                                                                                              ccgcAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                      >  3:1025063/1‑149 (MQ=255)
                                                                                              ccgcAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCaa                                                      <  1:809568/149‑1 (MQ=255)
                                                                                                gcAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGt                                                    >  3:353870/1‑149 (MQ=255)
                                                                                                          gcaATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACAt                                          <  2:1170430/149‑1 (MQ=255)
                                                                                                                 aTAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGAtt                                                                                                        >  1:284472/1‑80 (MQ=255)
                                                                                                                 aTAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACa                                           >  2:1209455/1‑141 (MQ=255)
                                                                                                                 aTAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACa                                           <  1:1209455/141‑1 (MQ=255)
                                                                                                                   aaCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGTTATTCCTCAAGTAAAAAAACATCTATTCCta                                 >  2:952097/1‑148 (MQ=255)
                                                                                                                   aaCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTg                                 >  3:248393/1‑149 (MQ=255)
                                                                                                                   aaCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTg                                 >  4:14501/1‑149 (MQ=255)
                                                                                                                    aCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAg                                                                                                                                                  >  3:44347/1‑35 (MQ=255)
                                                                                                                    aCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAg                                                                                                                                                  <  4:44347/35‑1 (MQ=255)
                                                                                                                    aCAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGc                                >  2:284223/1‑149 (MQ=255)
                                                                                                                      aGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTc                                    >  2:167438/1‑143 (MQ=255)
                                                                                                                      aGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGa                              >  2:649274/1‑149 (MQ=255)
                                                                                                                       gAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGAt                             <  1:1555606/149‑1 (MQ=255)
                                                                                                                         aaCCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGAttt                           >  1:213061/1‑149 (MQ=255)
                                                                                                                          aCCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGCTTtc                          >  4:1415215/1‑149 (MQ=255)
                                                                                                                          aCCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTtc                          <  2:213061/149‑1 (MQ=255)
                                                                                                                           ccTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTtct                         <  4:919980/149‑1 (MQ=255)
                                                                                                                             tGAAACAC‑AAAACGGCAGCCCTTGGGCTGCCGATTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTca                       >  4:313098/1‑149 (MQ=255)
                                                                                                                             tGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTca                       <  2:832323/149‑1 (MQ=255)
                                                                                                                                 acac‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaa                   <  1:649274/149‑1 (MQ=255)
                                                                                                                                   ac‑aAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaa                 >  3:202115/1‑149 (MQ=255)
                                                                                                                                    c‑aAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa                <  4:1164457/149‑1 (MQ=255)
                                                                                                                                    c‑aAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACaaaaaa                <  2:1074592/149‑1 (MQ=255)
                                                                                                                                      aaaaCGGCAGCTCAAGGGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACATAAAAg               >  3:897753/1‑149 (MQ=255)
                                                                                                                                      aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAg               >  2:1409536/1‑149 (MQ=255)
                                                                                                                                      aaaaCGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAg               >  1:872705/1‑149 (MQ=255)
                                                                                                                                          cGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTc           <  1:1569318/149‑1 (MQ=255)
                                                                                                                                            gCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGt         >  4:389243/1‑149 (MQ=255)
                                                                                                                                            gCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGt         >  4:1173232/1‑149 (MQ=255)
                                                                                                                                              aGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTg       >  2:956449/1‑149 (MQ=255)
                                                                                                                                                 ccTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGAt             >  3:60656/1‑140 (MQ=255)
                                                                                                                                                 ccTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGAt             <  4:60656/140‑1 (MQ=255)
                                                                                                                                                 ccTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACa    >  3:62410/1‑149 (MQ=255)
                                                                                                                                                 ccTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACa    >  1:191737/1‑149 (MQ=255)
                                                                                                                                                  cTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACaa   >  2:1007417/1‑149 (MQ=255)
                                                                                                                                                   ttGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACaaa  >  4:707561/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCCAGGCTACGCGTGATGCGCTGCCGCTGATGGAAAATGCATACCGTCAGATGCAGATGAATGCGCAAGCTGAAAAAGTAGCGAAAATCATCGCCGCAAACAGCAGCAATACATAACAGAAACCTGAAACAC‑AAAACGGCAGCCCTTGAGCTGCCGTTTTTTTATTCTGTCAGTTGTGAAACTGAAGCGATTTAGTCGCTATCGATCTCATCAAATATGGCTCGCTTTGAGATATTCCTCAAGTAAAAAAACATCTCTTCCTGCGATTTCTCACAAAAAAGATTCGTTGACAAA  >  NC_000913/2736767‑2737061

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑