Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,900,410 A→T 7.2% P56P (CCT→CCA yieL ← putative xylanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,900,4100AT 7.2% 756.6 / 37.2 221P56P (CCT→CCAyieLputative xylanase
Reads supporting (aligned to +/- strand):  ref base A (105/100);  new base T (11/5);  total (116/105)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-02

Only 100 of 226 total aligned reads displayed.

GTACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAG  >  NC_000913/3900262‑3900558
                                                                                                                                                    |                                                                                                                                                    
gTACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCa                                                                                                                                                      >  3:622209/1‑149 (MQ=255)
 tACACCATTAACATCGAAAAAATACTCGTACTGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGATTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAg                                                                                                                                                     <  3:1044570/149‑1 (MQ=255)
  acacCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAgg                                                                                                                                                    >  2:1050890/1‑149 (MQ=255)
  acacCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAgg                                                                                                                                                    <  2:374010/149‑1 (MQ=255)
  acacCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAgg                                                                                                                                                    <  4:1465784/149‑1 (MQ=255)
      cATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAg                                                                                                                                                     >  3:686422/1‑144 (MQ=255)
       aTCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaa                                                                                                                                               >  1:1319947/1‑149 (MQ=255)
        tCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaaa                                                                                                                                              <  2:1640371/149‑1 (MQ=255)
         cAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCTGGGGCaaaa                                                                                                                                             <  2:1503493/149‑2 (MQ=255)
         cAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAg                                                                                                                                             <  1:1644071/149‑1 (MQ=255)
             aTTGAAAAAATACTCGTACAGGTTGCCTTTCATGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCATGGGCAAAGTAGc                                                                                                                                         >  3:811155/1‑149 (MQ=255)
              ttGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCg                                                                                                                                        >  1:526888/1‑149 (MQ=255)
              ttGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCg                                                                                                                                        <  1:293143/149‑1 (MQ=255)
                  aaaaaTACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAg                                                                                                                                    <  4:1328017/149‑1 (MQ=255)
                      aTACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGTGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCa                                                                                                                                <  2:1275258/149‑1 (MQ=255)
                         cTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGc                                                                                                                             <  1:951837/149‑1 (MQ=255)
                         cTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGGCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGc                                                                                                                             >  4:881765/1‑149 (MQ=255)
                          tCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCt                                                                                                                            >  4:1266832/1‑149 (MQ=255)
                          tCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCt                                                                                                                            >  2:234014/1‑149 (MQ=255)
                            gTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACTCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGt                                                                                                                          <  3:186919/149‑1 (MQ=255)
                            gTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGt                                                                                                                          >  4:1619882/1‑149 (MQ=255)
                             tACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGtt                                                                                                                         >  2:1029415/1‑149 (MQ=255)
                                agtgtGCCTTTCAGGATAGGGGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGACCCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATc                                                                                                                      >  2:1446223/5‑149 (MQ=255)
                                aGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATc                                                                                                                      >  4:602296/1‑149 (MQ=255)
                                  gTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGt                                                                                                                                  <  3:1613252/135‑1 (MQ=255)
                                  gTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCgg                                                                                                                    >  4:1142026/1‑149 (MQ=255)
                                   ttGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTCCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGc                                                                                                                   >  4:724770/1‑149 (MQ=255)
                                      ccTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTCGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCGCCACTACAAAGACATTTGTTGCCCCAGGGGCAAAGTAGCGAAAGGTCTCGCTGTTATCGGCAtt                                                                                                                >  4:1007754/1‑149 (MQ=255)
                                        tttCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGa                                                                                                              >  3:1184877/1‑149 (MQ=255)
                                         ttCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGAc                                                                                                             <  3:1385067/149‑1 (MQ=255)
                                         ttCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGAc                                                                                                             <  3:230337/149‑1 (MQ=255)
                                           cAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCt                                                                                                           <  2:749894/149‑1 (MQ=255)
                                             ggATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGa                                                                                                         <  3:425282/149‑1 (MQ=255)
                                             ggATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGa                                                                                                         <  1:89136/149‑1 (MQ=255)
                                             ggATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGa                                                                                                         <  1:11240/149‑1 (MQ=255)
                                              gATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGAc                                                                                                             >  4:609917/1‑144 (MQ=255)
                                              gATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAg                                                                                                        >  1:1388253/1‑149 (MQ=255)
                                               aTAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGt                                                                                                       >  3:483959/1‑149 (MQ=255)
                                               aTAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGt                                                                                                       <  1:910429/149‑1 (MQ=255)
                                                   gtgtGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGAAGCGaaa                                                                                                                                     >  4:16765/1‑115 (MQ=255)
                                                         gcCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGt                                                                                             <  3:286865/149‑1 (MQ=255)
                                                              gACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCAc                                                                                        >  3:642291/1‑149 (MQ=255)
                                                                 cAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCATCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGATTCACATACTGTTTCACCgg                                                                                     <  2:1118735/149‑1 (MQ=255)
                                                                      tCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaaa                                                                                                                                              >  2:965261/1‑87 (MQ=255)
                                                                      tCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCAtt                                                                                                                <  2:940415/117‑1 (MQ=255)
                                                                      tCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCAtt                                                                                                                >  1:940415/1‑117 (MQ=255)
                                                                      tCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAg                                                                                >  3:1063173/1‑149 (MQ=255)
                                                                      tCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAg                                                                                >  3:1021084/1‑149 (MQ=255)
                                                                      tCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAg                                                                                >  2:1342324/1‑149 (MQ=255)
                                                                       ccTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGt                                                                               <  1:1615085/149‑1 (MQ=255)
                                                                       ccTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGt                                                                               <  1:1442499/149‑1 (MQ=255)
                                                                        cTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATa                                                                                                 >  1:767687/1‑130 (MQ=255)
                                                                        cTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTg                                                                              >  1:250583/1‑149 (MQ=255)
                                                                        cTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTg                                                                              <  1:1275070/149‑1 (MQ=255)
                                                                        cTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTg                                                                              >  2:311410/1‑149 (MQ=255)
                                                                         tGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTgg                                                                             >  4:660976/1‑149 (MQ=255)
                                                                                tCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCgg                                                                      >  1:1085486/1‑149 (MQ=255)
                                                                                 cTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCg                                                                                      <  2:1387167/132‑1 (MQ=255)
                                                                                  tttGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaaa                                                                                                                                              <  3:35071/75‑1 (MQ=255)
                                                                                  tttGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCaaa                                                                                                                                              >  4:35071/1‑75 (MQ=255)
                                                                                  tttGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggc                                                                    >  3:86849/1‑149 (MQ=255)
                                                                                       tCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggctggct                                                               >  2:284477/1‑149 (MQ=255)
                                                                                        cATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGgctggctg                                                              <  4:229218/149‑1 (MQ=255)
                                                                                         aTCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg                                                             <  4:235016/149‑1 (MQ=255)
                                                                                         aTCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGctggctgg                                                             <  1:866255/149‑1 (MQ=255)
                                                                                          tCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCtggctggc                                                            >  1:616729/1‑149 (MQ=255)
                                                                                            gggTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCat                                                          >  4:190853/1‑149 (MQ=255)
                                                                                              gTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCatat                                                        >  2:1077472/1‑149 (MQ=255)
                                                                                               tGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATc                                                       <  2:1602356/149‑1 (MQ=255)
                                                                                                   tatTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggc                                                                    <  4:516567/132‑1 (MQ=255)
                                                                                                   tatTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGggc                                                                    >  3:516567/1‑132 (MQ=255)
                                                                                                    atTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCgctg                                                  <  2:551955/149‑1 (MQ=255)
                                                                                                           ggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGc                                           >  3:1223857/1‑149 (MQ=255)
                                                                                                           ggaaccggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGc                                           >  2:1548853/1‑149 (MQ=255)
                                                                                                              accggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCa                                        <  1:1548853/149‑1 (MQ=255)
                                                                                                                cggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACt                                      >  1:324847/1‑149 (MQ=255)
                                                                                                                  gaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGa                                    <  4:483959/149‑1 (MQ=255)
                                                                                                                        ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCa                                                                                                                                >  2:774641/1‑51 (MQ=255)
                                                                                                                        ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCAcc                                                                                       >  4:866204/1‑92 (MQ=255)
                                                                                                                        ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGaa                              >  1:150241/1‑149 (MQ=255)
                                                                                                                        ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGaa                              >  1:1055246/1‑149 (MQ=255)
                                                                                                                        ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGaa                              >  4:1202335/1‑149 (MQ=255)
                                                                                                                        ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGaa                              >  1:98840/1‑149 (MQ=255)
                                                                                                                        ccaccaCTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGCCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGaa                              >  3:1132813/1‑149 (MQ=255)
                                                                                                                         caccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGaaa                             <  1:1029415/149‑1 (MQ=255)
                                                                                                                          accacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAAt                            >  3:1537349/1‑149 (MQ=255)
                                                                                                                               tACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGc                       >  1:681645/1‑149 (MQ=255)
                                                                                                                                  agagACATTTTTTGCCACAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGa                                    <  4:630793/133‑1 (MQ=255)
                                                                                                                                    agaCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCa                  <  2:865156/149‑1 (MQ=255)
                                                                                                                                       cATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATa               <  3:948663/149‑1 (MQ=255)
                                                                                                                                        ctactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAg              >  4:398992/5‑149 (MQ=255)
                                                                                                                                         tactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGaaa                             >  4:1164714/4‑133 (MQ=255)
                                                                                                                                         tactttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGc             >  3:1308110/4‑149 (MQ=255)
                                                                                                                                          actttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGcc            >  3:258498/3‑149 (MQ=255)
                                                                                                                                           ttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCa           <  3:602296/149‑1 (MQ=255)
                                                                                                                                            tttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCgg                                                                      <  3:1669408/89‑1 (MQ=255)
                                                                                                                                            tttGCCCCTGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAg          >  1:855543/1‑149 (MQ=255)
                                                                                                                                                  ccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaa    <  4:1537349/149‑1 (MQ=255)
                                                                                                                                                  ccAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaa    >  2:949873/1‑149 (MQ=255)
                                                                                                                                                    aGGGGCAAAGTAGCGAAAGGTAACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATACGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGCTaacg  >  1:226202/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTACACCATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAG  >  NC_000913/3900262‑3900558

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