Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,000,064 T→A 6.0% I19F (ATC→TTC)  PP_0861 ← outer membrane ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,000,0640TA 6.0% 199.6 / 2.4 66I19F (ATC→TTC) PP_0861outer membrane ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (39/23);  new base A (2/2);  total (41/25)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CTTCTGCGAGGCGGCCTTCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTACCATTCGTGAA  >  NC_002947/999921‑1000206
                                                                                                                                                 |                                                                                                                                                 
cTTCTGCGAGGCGGCCTTCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGTCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATggcag                                                                                                                                              >  3:253143/1‑149 (MQ=255)
 ttCTGCGAGGCGGCCTTCTCCACCTTGTAATCGGTGCCGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATggcagg                                                                                                                                             >  4:397311/1‑149 (MQ=255)
   cTGCGAGGCGGCCTTCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATggcaggc                                                                                                                                            <  3:398409/148‑1 (MQ=255)
      cGAGGCGGCCTTCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGa                                                                                                                                        <  5:141216/149‑1 (MQ=255)
      cGAGGCGGCCTTCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGa                                                                                                                                        >  6:62292/1‑149 (MQ=255)
          gcggcCTTCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAg                                                                                                                                    >  1:94418/1‑149 (MQ=255)
            ggcCTTCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGc                                                                                                                                  <  5:19478/149‑1 (MQ=255)
                ttCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTgg                                                                                                                              >  6:235273/1‑149 (MQ=255)
                     cACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGaa                                                                                                                       <  5:62292/148‑1 (MQ=255)
                       ccTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCa                                                                                                                                     <  4:251798/135‑1 (MQ=255)
                       ccTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCa                                                                                                                                     >  3:251798/1‑135 (MQ=255)
                         ttGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTca                                                                                                                   >  3:316684/1‑148 (MQ=255)
                         ttGTAATCGGTGATGGCTTGTTCAGCCTTTAAGTTGACGCTGGTGGCATCGAGAGTGATTGCGCTTTTGGAGGACGGAAAGTCGGCATATACAGACTAGGCACTCAGGGCGACAACGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAAATca                                                                                                                   >  7:191085/1‑148 (MQ=255)
                          tGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTcaca                                                                                                                 >  7:236391/1‑149 (MQ=255)
                              aTCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGCTCAGGCGTGG‑‑‑CGAACTCAcag                                                                                                                >  6:5248/1‑146 (MQ=255)
                               tccgTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCGTATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCAcca                                                                                                                                        >  4:222405/4‑122 (MQ=255)
                               tCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCGTATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGa                                                                                                                                        <  3:222405/124‑1 (MQ=255)
                                 ggggCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGAt                                                                                                           <  4:253143/145‑1 (MQ=255)
                                 ggTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGAt                                                                                                           <  8:236391/148‑1 (MQ=255)
                                  gTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATgg                                                                                                         >  2:192699/1‑149 (MQ=255)
                                     cTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGgtac                                                                                                      <  2:94418/149‑1 (MQ=255)
                                        gCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGgtacgta                                                                                                   >  1:367702/1‑149 (MQ=255)
                                            gTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTg                                                                                                <  7:387644/148‑1 (MQ=255)
                                             ttCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGgc                                                                                              >  6:140194/1‑149 (MQ=255)
                                                aGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGaa                                                                                                                       <  1:246748/122‑1 (MQ=255)
                                                aGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGaa                                                                                                                       >  2:246748/1‑121 (MQ=255)
                                                    ttGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACAGATCGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTAACAGCAGATGGTACGTACTGGCGCATTg                                                                                        >  1:131064/1‑149 (MQ=255)
                                                     tGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGc                                                                                       >  1:30364/1‑148 (MQ=255)
                                                       ccGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTca                                                                                                                   <  4:363248/118‑1 (MQ=255)
                                                       ccGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTca                                                                                                                   >  3:363248/1‑118 (MQ=255)
                                                          ttGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGaa                                                                                 <  7:92487/149‑1 (MQ=255)
                                                            gACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAAc                                                                                <  2:367702/148‑1 (MQ=255)
                                                                        aTCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTa                                                                   >  4:314170/1‑149 (MQ=255)
                                                                        aTCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTa                                                                   >  8:243932/1‑149 (MQ=255)
                                                                        aTCGAGAGTGATGGCGATGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGtt                                                                   >  8:59570/1‑148 (MQ=255)
                                                                         tCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTAt                                                                  >  6:374646/1‑149 (MQ=255)
                                                                            agagTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTAtt                                                                 >  8:213730/1‑147 (MQ=255)
                                                                             gagTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAg                                                              >  5:25667/1‑149 (MQ=255)
                                                                                tGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCg                                                           >  6:43691/1‑149 (MQ=255)
                                                                                tGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGAAGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCg                                                           >  1:122941/1‑149 (MQ=255)
                                                                                 gATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCgg                                                          <  4:316684/149‑1 (MQ=255)
                                                                                  aTGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGt                                                         <  1:192699/149‑1 (MQ=255)
                                                                                    ggCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTg                                                        >  2:13375/1‑148 (MQ=255)
                                                                                        cTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCt                                                   <  3:314170/149‑1 (MQ=255)
                                                                                         tGTTGGCGGCCTGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACCGGCGCATTGCAAAGCACCTTCCCCAAGGTATTGAGGCGGTGGAGCt                                                   >  6:97417/1‑148 (MQ=255)
                                                                                             ggcggcCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACAGCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACCGGCGCATTGCAAAGAACCTTCCCAAAGGTATTGAGGCGGTGGATCTTAAt                                               <  5:97417/148‑1 (MQ=255)
                                                                                              gcggcCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGAAGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGt                                             <  2:122941/149‑1 (MQ=255)
                                                                                              gcggcCGGATCGGCGCCATATACCGACGAGGCACTCAGCGCGACACCTATGGCAGCCTCCATCACGCGTGGCTCCAAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTa                                                 <  1:354374/148‑1 (MQ=255)
                                                                                                            ggcATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGAc                                >  5:262156/1‑148 (MQ=255)
                                                                                                               atatACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTaa                             >  7:15447/1‑148 (MQ=255)
                                                                                                                     cgacgaGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAAt                                               <  8:145445/124‑1 (MQ=255)
                                                                                                                     cgacgaGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAAt                                               >  7:145445/1‑124 (MQ=255)
                                                                                                                     cgacgaGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAAtt                        >  3:124055/1‑147 (MQ=255)
                                                                                                                     cgacgaGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAAtt                        >  3:214420/1‑147 (MQ=255)
                                                                                                                      gacgaGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAAtt                        >  4:125906/1‑146 (MQ=255)
                                                                                                                       acgaGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTga                      >  5:280914/1‑147 (MQ=255)
                                                                                                                               aCTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTa             <  6:25667/149‑1 (MQ=255)
                                                                                                                                cTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTAc            >  2:12436/1‑149 (MQ=255)
                                                                                                                                   aGGGCGACACCGAAGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTgaga                    <  6:124583/138‑1 (MQ=255)
                                                                                                                                   aGGGCGACACCGAAGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTgaga                    >  5:124583/1‑138 (MQ=255)
                                                                                                                                     ggCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTACCATTc       >  1:167495/1‑149 (MQ=255)
                                                                                                                                     ggCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTACCATTc       <  5:140194/149‑1 (MQ=255)
                                                                                                                                     ggCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTACCATTc       >  8:238446/1‑149 (MQ=255)
                                                                                                                                        gACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTACCATTCGt     >  1:165324/1‑148 (MQ=255)
                                                                                                                                          cacCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTACCATTCGTg    >  3:84074/1‑147 (MQ=255)
                                                                                                                                          cacCGATGGCAGGCACGATAAGGCGTGG‑‑‑CTAACTCACAGCAGTTGTTACGTACTGGCTCATTGCAAAGAACCTTCCCCAAGGTTTTTAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACATTTGAGAGTCGTTACCATTCGTgaa  <  7:59570/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTTCTGCGAGGCGGCCTTCTCCACCTTGTAATCGGTGCTGGCTTGTTCAGCCTTGCCGTTGACGCTGGTGGCATCGAGAGTGATGGCGCTGTTGGCGGCCGGATCGGCGGCATATACCGACGAGGCACTCAGGGCGACACCGATGGCAGGCACGATCAGGCGTGG‑‑‑CGAACTCACAGCAGATGGTACGTACTGGCGCATTGCAAAGAACCTTCCCCAAGGTATTGAGGCGGTGGATCTTAATGTAAGCAATTTTGACTAACAATTGAGAGTCGTTACCATTCGTGAA  >  NC_002947/999921‑1000206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: